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1 INTRODUCTION 



In the Standard Model (SM), CP violation arises from a single complex phase in the Cabibbo— 
Kobayashi-Maskawa (CKM) quark-mixing matrix [1]. Decays of B mesons into charmless hadronic 
final states with three kaons are dominated by 6 ^ sss gluonic penguin amplitudes, with smaller 
contributions from electroweak penguins, while other SM amplitudes are suppressed by CKM fac- 
tors [2]. The time-dependent CP-asymmetry is obtained by measuring the proper time difference 
At = tcp — ^tag between a fully reconstructed neutral B meson {Bqp) in the final state K~ K^^ 
and the partially reconstructed recoil B meson (-Btag)- decay products of -B^ag Provide ev- 
idence that it decayed either as B^ or B^ (flavor tag). The decay rate f+(f-) when the tagging 
meson is a is given by 

„-|At|/TaO 

f±(At) = — [ l±5sin(AmdAt)^ Ccos(AmdAt)], (1) 

where t^o is the neutral B meson mean lifetime and Amd is the B^-B*^ oscillation frequency. The 
parameters C and S describe the magnitude of CP violation in the decay and in the interference 
between decay and mixing, respectively. The time-dependent CP-violating asymmetry is defined as 
Acp = (f+ — f_)/(f+-|-f_). In the SM, we expect C = because there is only one decay mechanism 
and direct CP violation requires amplitudes with different phases. Neglecting CKM-suppressed 
contributions and assuming that K^K~K^ decay proceeds through an S-wave, leading to a CP- 
odd final state, the time-dependent CP-violating parameter S in this decay and B^ J/t/j are 
both equal to the same parameter sin 2/3 [3], where the latter decay is dominated by tree diagrams. 
Since many scenarios of physics beyond the SM introduce additional diagrams with heavy particles 
in the penguin loops and corresponding new phases, comparison of CP- violating observables with 
SM expectations is a sensitive probe for new physics [4]. Measurements of sin2/3 in B decays to 
charmonium such as B^ J/ipKg have been reported by the BABAR [5] and Belle [6] collaborations, 
and the world average for sin2/3 (0.736 ± 0.049 [7]) is in good agreement with SM expectations [8]. 
A deviation from this value in the case of loop-dominated channels might signal the presence of 
physics beyond the SM. 

Measurements of the CP asymmetry in the decays B^ (t'K^ and B^ currently 
have large statistical uncertainties [9,10]. More accurate CP asymmetry measurements have been 
performed in the final state K'^K^Kg, (excluding B^ — > (f'^s) [H]) which has a branching fraction 
several times larger than the resonant modes, The CP content of the final state, which is a priori 
unknown, is estimated using an angular-moment analysis. 

In this document we report preliminary measurements of the time dependent CP asymmetry 
in the CP conjugate state B^ K+R-R^. 

2 THE BABAR DETECTOR AND DATASET 

This measurement is based on a sample of approximately 227 million BB pairs collected at the 
T{4S) resonance with the BABMl detector [12] at the PEP-II asymmetric-energy e'^e~ storage 
ring [13] located at the Stanford Linear Accelerator Center. The BABAR detector is fully described 
elsewhere [12]. The detector systems used in this analysis are a charged-particle tracking sys- 
tem consisting of a five-layer silicon vertex tracker (SVT) and a 40-layer drift chamber (DCH) 
surrounded by a 1.5-T solenoidal magnet with an instrumented flux return (IFR), an electromag- 
netic calorimeter (EMC) composed of 6580 CsI(Tl) crystals, and a detector of internally reflected 
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Cherenkov light (DIRC) providing excellent charged K — n identification up to a momentum of 
4.5 GeV/c, which is the relevant momentum range for this analysis. 

3 EVENT RECONSTRUCTION 

The K~ candidate (Bcp) is reconstructed by combining a pair of oppositely charged 

tracks extrapolated from a common vertex and a candidate. For the charged tracks we require 
at least 12 measured drift-chamber coordinates and a minimum transverse momentum of 0.1 GeV/c. 
The tracks must also originate within ±10 cm along the beam axis and 1.5 cm in the transverse 
plane, with respect to the nominal beam spot. Charged kaons are distinguished from pion and 
proton tracks via a requirement on a likelihood ratio that combines dE/dx information from the 
SVT and the DCH for tracks with momentum p < 0.7GeV/c. For tracks with higher p, dE/dx in 
the DCH and the Cherenkov angle and the number of photons as measured by the DIRC are used 
in the likelihood. These particle identification criteria limit the rate of pion misidentification as a 
kaon to less than 2%, with an efficiency of 70%. 

We identify a candidate as in the BABAR analysis of the decay B^ J/ipK^ analysis [14] 
either as a cluster of energy deposited in the EMC or as a cluster of hits in two or more layers of the 
IFR that cannot be associated with any charged track in the event. The energy is not measured, 
therefore, we determine the laboratory momentum from its flight direction as measured from 
the EMC or IFR cluster, and the constraint that the invariant K~ mass agree with the 
known B^ mass. In those cases where the is detected in both the IFR and EMC we use the 
angular information from the EMC, because it has higher precision. In order to reduce background 
from vr'^ decays, we reject an EMC candidate cluster if it forms an invariant mass between 100 
and 150 MeV/c^ with any other neutral cluster in the event under the 77 hypothesis, or if it has 
energy greater than 1 GeV and contains two shower maxima consistent with two photons from a 
vr'^ decay. The remaining background of candidates due to photons and overlapping showers 
is further reduced with the use of a neural network MM emc- The MM emc is constructed from 
cluster shape variables, trained using as signal measured B^ — > J/ipK^ events and as background 
measured K~ events which lie outside the signal region, and tested on e+e~ — > KgK^)^ 
events. 

The results are extracted using an extended unbinned maximum likelihood fit. We parameterize 
the distributions of kinematic and topological variables for signal and background events in terms 
of probability density functions (PDFs) [17]. The selection requirements for these variables is loose 
to include background dominated regions which can then be extrapolated into the signal region. 
The main source of background, estimated from data, comes from random combinations of tracks 
produced in events of the type e+e~ — > qq, where q = u,d,s,c (continuum). Background from 
decays of B mesons in other final states with and without charm is estimated using Monte Carlo 
simulation. 

In the following we describe the event variables used in the maximum likelihood fit to character- 
ize the signal and background: the energy difference AE = E^ — ^^/s, where E^ is the energy of the 
B candidate and y/s is the total energy, both evaluated in the T(4S') rest frame, a neural network 
MM built with topological quantities, and At, described in Section ^ For signal events, AE is 
expected to peak at zero, with a broad tail for positive values of AE. We require AE < 0.08 GeV, 
in order to be able to fix the shape of background under the signal peak. The AE resolution is 3.0 
MeV, which includes the different resolutions for EMC and IFR events. This resolution has been 
validated on data using reconstructed B^ J/t/jK^ events. 
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Continuum events are characterized by a jet-like topology in the T(45') rest frame, because 
of the large amount of phase space in the decay, while B mesons are produced almost at rest, 
so particles produced in B decays are distributed isotropically. One can then define a set of 
topological variables to quantify the sphericity of the event. One such quantity is the angle OsPH 
between the sphericity axis of the Bcp candidate and the sphericity axis formed from the other 
charged and neutral particles in the event. We also use the cosine of the angle 6b between the Bcp 
momentum and the beam axis, and the sum of the momenta pi of the other charged and neutral 
particles in the event weighted by the Legendrc polynomials Lo(0i) and L2(^i) where 0i is the angle 
between the momentum of particle i and the thrust axis of the Bcp candidate. We use other 
variables characterizing a final state with a K^. One is the reconstructed energy difference AE^is, 
calculated in the laboratory frame as the difference between total reconstructed energy of the event, 
where the energy of the neutral particles is calibrated for an electromagnetic shower, and the two 
reconstructed kaon energies. The other is the cosine of the polar angle of the missing momentum, 
Pmiss, calculated as the difference between the sum of beam momenta and all tracks and EMC 
clusters, in the laboratory frame, excluding the candidate. We combine these variables in a 
neural network A^A^, which peaks at for continuum events and at 1 for signal events. We apply 
a selection that retains 80% of signal events and rejects 84% of continuum background. The rest 
of the events are used for the maximum likelihood fit. 

The remaining background originates from B decays where a neutral or charged pion is missed 
during the reconstruction (peaking B background). Since the branching fractions for these decay 
modes (5° ^ K+ R- K*^{Kl7r'^), B^ K*+ {K+tt^)K^ Kl, B+ K+ R- K*+{KIt:+), and 
B^ — > K^K*~ (^K~ TT^)) are not known, we build a cocktail of exclusive Monte Carlo samples, 
weighted with the relative efficiency, and the yield is floated in the final fit. The rest of the 
background originating from B decays comes from the combinations of particles originating from 
both B mesons that have continuum-like values of /\E, so these events are included in the continuum 
component by the fit, without generating a bias in the fitted values of S and C parameters^. 

We suppress background from B decays that proceed through ah ^ c transition leading to the 
K~^K~K^ final state by applying invariant mass cuts to remove J/0, Xdd-> and ip{2S) decaying 
into K~^K~, and and Df decays into K^K^. 

All the other tracks and clusters that are not associated with the reconstructed B^ K~ 
decay are used to form the -B^agi flavor is determined with a multivariate tagging algorithm [15]. 
The tagging efficiency e and mistag probability w in five hierarchical and mutually exclusive cat- 
egories are measured using fully reconstructed B'^ decays into the D(*)^X+ {X~^ = vr"*", p"*", a^) 
and J/ipK*^ {K*^ K~^Tr^) fiavor eigenstates (-Bg^y sample). The analyzing power e(l — 2w)'^ is 
(30.3 ± 0.4)%. 

A detailed description of the At reconstruction algorithm is given in Ref. [14]. 

4 MAXIMUM LIKELIHOOD FIT 

The CP asymmetry parameters are extracted from a K'^K~Kl sample which excludes B° (l)K^ 
with an invariant mass veto: \m{K'^ K~) — m{(f))\ > 15MeV/c^. This excludes B^ — > (fjK^ events by 
three standard deviations. The average Az resolution is 190/im, dominated by the tagging vertex in 
the event. Thus, we can characterize the resolution using the much larger Bq^^ sample, which we 
use for signal parameterization. The amplitudes for the Bcp asymmetries and for the Bq^^ fiavor 

''This conclusion has been obtained simulating the fit on a sample of fully simulated Monte Carlo events. 
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oscillations are reduced by the same factor due to wrong tags. Both distributions are convolved with 
a common At resolution function. Backgrounds are accounted for by adding terms to the likelihood, 
incorporated with different assumptions about their At evolution and resolution function [14]. The 
At resolution function is parameterized as a sum of two Gaussian distributions with different mean 
values whose widths are given by a scale factor times the event-by-event uncertainty a At- A third 
Gaussian distribution, with a fixed large width, accounts for a small fraction of outlying events [5]. 
For the time-dependent fit we retain events that have |At| < 20 ps and whose estimated uncertainty 
(TAt is less than 2.5 ps. 

Since we measure the correlation among the observables to be small in the data sample used in 
the fit (the largest is 2.9 % between AE and MM) we take the probability density function Vj^ for 
each event j to be a product of the PDFs for the separate observables. For each event hypothesis 
i (signal, backgrounds) and tagging category c, we define Vj^ = Vi{AE) ■ Vi{MM) ■ Vi{At; a^t, c). 
The likelihood function for each decay chain is then 



where Ni^c is the yield of events of hypothesis i obtained from the fit in category c, and Nc is the 
number of category c events in the sample. The total sample consists of 77577 K~^K~K^ candidates. 
The total reconstruction efficiency is (e) = (23.1 ± 0.6)%. We fixed in the fit SB-bkg = 0.42 
(as estimated from full Monte Carlo simulation of generic neutral and charged S-decays) and 
Cs-bkg = (as the SM expectations). From the fit we find 777 ± 80 K^K~K^ signal events 
(B^ — > (j)K^ excluded). The signal yield agrees with the branching fraction determined in the 
K+K'Kl final state within one standard deviation, but the uncertainty in efficiency is large. 
Figure^shows the AE distribution together with the result from the fit after a cos 9h cut to enhance 
signal, where Oh is the angle between the candidate and the parent Bqp flight direction in 
the K~^K~ rest frame. The neural network output distribution together with the result from the 
fit after a requirement on the likelihood to enhance the sensitivity. The fit was tested with both a 
parameterized simulation of a large number of data-sized experiments and a full detector simulation. 
The likelihood of our data fit agrees with the likelihoods from fits to the simulated data. Figure |2l 
shows the comparison between data and Monte Carlo simulated events of the signal to background 
likelihood ratio -^sig/(-^^sig + -^bkg) distribution. It shows the goodness of the agreement between 
data and Monte Carlo parameterization event-by-event. The fit was also verified with our J/ipK^ 
data sample to check the fitted central value of sin2/? and the sign of the CP eigenstate definition. 



The measurement of the CP content has been done in the K~^K~Kg final state from an angular 
moments analysis [11]. Since the K^K^K^ final state has higher purity than can be obtained in 
the final state, the angular analysis has not been repeated with our sample, but the results on 
the CP-conjugate state have been used. In order to take into account differences in the efficiency 
across the Dalitz plot between the K^K^K^ and K^K^K^ samples, which can change the relative 
amount of CP-even(odd) fraction in different m{K^K^) regions, we use the feven fraction measured 
in seven m{K~^ K~) bins (excluding the (j) region) in B° K+K-Kl and compute a re-weighted 
average using yields in the same mass bins. Tabled shows the K~ efficiencies, 

yields and the measured K^K^K^ feven- Variations of the signal yield are shown also in the upper 
plot of Fig. H 





5 ESTIMATION OF CP CONTENT 
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Figure 1: Distribution of the event variable /S.E after a cos Oh cut (left) and neural network output 
after a requirement on the likelihood, calculated without the plotted variable. The signal efficiency 
for the selection and likelihood requirements is 31% for /S.E and 8% for the neural network output. 
The solid line represents the fit result for the total event yield and the dashed line for the total 
background. The dotted line represents the continuum background, only. 



Table 1: Average efficiencies, yields for JB° — K^K K^, and fevenj calculated in the K^K 
final state, for seven K^K~ mass bins (excluding ^ region). 



m{K+K- 


-) GeV/c2 


{e) 


Signal yield 




[1.1 


1.3] 


0.153 


67.6 ± 20.8 


1.10 ± 0.18 


[1.3 


1.5] 


0.172 


44.0 ± 21.3 


0.99 ± 0.14 


[1.5 


1.9] 


0.188 


93.9 ± 28.0 


0.93 ± 0.21 


[1.9 


2.3] 


0.223 


146.4 ± 28.1 


0.95 ± 0.16 


[2.3 


2.7] 


0.258 


117.4 ± 24.0 


0.79 ± 0.24 


[2.7 


3.1] 


0.271 


93.7 ± 21.8 


0.74 ± 0.25 


[3.1 


4.9] 


0.242 


141.9 ± 37.2 


0.96 ± 0.43 



Out of the (j) region the sample consists mainly of S-wave decays, giving an fraction close 
to 1. We find the total fraction of CP-odd final states: 

fodd = 0.92 ± 0.07 (stat) ± 0.06 (syst), 

where the systematic uncertainty is evaluated in the BABMi — > K^K~K^ analysis [11]. 

6 SYSTEMATIC STUDIES 

We consider systematic uncertainties in the CP coefficients S and C due to the parameterization of 
PDFs for the event yield in signal and background by varying the parameters within one standard 
deviation (evaluated from a fit to Monte Carlo simulated events). Since the real CP content of the B 
background is not known, we vary SB-bkg and CB-bkg in ^ conservative interval. About 50% of these 
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Figure 2: Distribution of signal to background likelihood ratio -^gig/l^gig + -^bkg)' '^^^ solid line 
represents the Monte Carlo simulation of the entire sample, and the dashed line the simulation of 
background only. 



events come from charged B decays, which show only direct CP-violation, while neutral B decays 
can violate CP both in mixing and decay. We therefore vary the B background CP parameters in the 
interval —0.5 < SB-bkgiCs-bkg) < 0.5, which corresponds to a uniform variation of the parameters 
for neutral B background in the whole physically allowed interval — 1 < SB-bkgiCs-bkg) < 1- We 
evaluate the uncertainty associated with the assumed parameterization of the At resolution function 
for signal and B background, a possible difference in the efficiency between reconstructed B^ and 
B^ decays, and the fixed values for Am^ and r^o, by varying the parameters within one standard 
deviation (extracted from a fit to the B^^^ sample). We also estimate uncertainties coming from 
possible SVT layer misalignments. The bias in the coefficients due to the fit procedure is included 
in the uncertainty without making corrections to the final results. Finally, we estimate the errors 
due to the effect of doubly CKM-suppressed decays [16]. We add these contributions in quadrature 
to obtain the total systematic uncertainty. The summary is reported in Table |21 



7 RESULTS 

The coefficients of the time-dependent CP asymmetry in B^ — > decays (excluding the 

B^ — > (t)K^ final state) determined by the fit are: 

S = 0.07 ± 0.28 (stat)l:|]:l^(syst), 
C = 0.54 ± 0.22 (stat)t[!;j]^,(syst). 

Figure shows the At distributions of the B^- and the i?'^-tagged subsets together with the raw 
asymmetry, with the result of the combined time-dependent CP-asymmetry fit superimposed. We 
also fit the CP parameters in the same m{K~^ K~) regions used to estimate the CP-odd fraction. 
The results are shown in Figure 01 The presence of both P- and S-wave decays in our CP sample 
dilutes the measurement of the sine coefficient. If we account for the measured CP-odd fraction, we 
can extract the SM parameter sin 2(3. Using the estimate of the CP content based on the angular 
moments, and setting C = in the fit, we get 

sin2/? = S/{2fodd - 1) = 0.09 ± 0.33 (stat)l|j:J^(syst) ± 0.10(syst CP-cont) 



13 




Figure 3: Distributions of At for —> K^K^K^ candidates with (top) B^- and (middle) i?'^-tags. 
The sohd hnes refer to the fit for all events; the dashed lines correspond to the background. The 
bottom plot shows the asymmetry. A requirement on signal-to-background ratio to enhance the 
signal is applied. 
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Figure 4: Distribution of signal yield (top), C (middle) and S (bottom) parameters in 7 different 
m{K~^K~) intervals. 
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Table 2: Summary of systematic uncertainties on the parameters S and C. The total systematic 
errors are obtained by adding in quadrature all individual sources. 



*M^l 1 YT'Ci 
OUUi Ct/ 




IS D{ ) 




A Ci—) 


Arud 


0.004 


-0.001 


0.000 


-0.001 


tbo 


0.01 


0.01 


0.00 


-0.00 


At model 


0.02 


0.02 


0.02 


0.01 


Tagging 


0.04 


0.04 


0.03 


0.03 


B background CP 


0.10 


0.10 


0.06 


0.07 


Signal and background PDFs 


0.03 


0.03 


0.02 


0.02 


fit biases 


0.00 


0.00 


0.02 


0.02 


SVT local alignment 


0.01 


0.01 


0.00 


0.00 


doubly-CKM-suppressed decays 


0.00 


0.00 


0.01 


0.01 


Total 


0.11 


0.12 


0.08 


0.09 



where the last error is due to uncertainty on the CP content. Since this uncertainty is multiplicative 
and the fitted value of S is close to 0, we conservatively computed this uncertainty shifting the 
measured value of S within 1 standard deviation. 

8 SUMMARY 

In a sample of 227 million BB mesons, we have obtained preliminary measurements of the CP 
content and CP parameters in the K^K^K^ final state that cxchides B^ 4>K^l decays. We 
estimated the fraction of S-wave events (CP-odd fraction) from measured value in the K^K^K^ 
final state, which has higher purity than K'^K'K^. The result shows the dominance of CP-odd 
final states. We compute the average of sin2/3g£jp and C parameter in K~^K~Kg and K~^K~ 
final states treating the systematic errors and the uncertainty on the CP content as completely 
correlated. This gives the most conservative estimation of the combined uncertainty. We obtain: 

[sm2p gff]^^ = 0.41 ± 0.18 (stat) ± 0.07 (syst) ± 0.11 (CPcontent) 
[C]av = 0-23 ± 0.12 (stat) ± 0.07 (syst) 

The result agrees within one standard deviation with the value of sin2/? in the 5° — (cc)if° 
decays [15]. 
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